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ABSTRACT 

We study the connections between the sun's convection zone and the evolution of the solar wind 
and corona. We let the magnetic fields generated by a 2.5D axisymmetric kinematic dynamo code 
(STELEM) evolve in a 2.5D axisymmetric coronal isothermal MHD code (DIP). The computations 
cover an 11 year activity cycle. The solar wind's asymptotic velocity varies in latitude and in time in 
good agreement with the available observations. The magnetic polarity reversal happens at different 
paces at different coronal heights. Overall sun's mass loss rate, momentum flux and magnetic braking 
torque vary considerably throughout the cycle. This cyclic modulation is determined by the latitudinal 
distribution of the sources of open flux and solar wind and the geometry of the Alfven surface. Wind 
sources and braking torque application zones also vary accordingly. 

Subject headings: solar dynamo — solar wind — Sun: corona — Sun: magnetic fields — Sun: con- 
vective zone 
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1. INTRODUCTION 

The sun's magnetic field varies in time following a 22- 
year cycle. The most visible manifestation of this peri- 
odic behaviour is the well-known 11 years sunspot cycle, 
the sunspots themselves corresponding to high local con- 
centrations (up to about 20 Mm wide) of vertical mag- 
netic field at the photosphere. The sun's photosphere is 
in reality permeated by magnetic fiux at several scales, 
and at all latitudes and longitudes. This fact brings for- 
ward the idea that the magnetic field plays a role in cou- 
pling the dynamics of the inner layers of the sun and its 
atmosphere. Studies about the convection zone (here- 
after CZ) and atmospheric magnetic fields have mostly 
been done separately, nonetheless. The first reason for 
this is that the photosphere essentially separates two re- 
gions of plasma with different regimes (/3 ^ 1 and /3 ^ 1) 
and disparate characteristic time and length scales, mak- 
ing any numerical investigations comprising them both a 
real challenge. The second reason is that the actual cou- 
pling and/or cross-transport phenomena are still poorly 
understood (except perhaps at very small-scale, or in 
simple scenarios). Finally, the third reason is linked to 
the use of different observational techniques for gathering 
data relating to the dense layers below the optically thick 
photosphere and to the essential ly optically thin a nd rar- 
efied atmosphere (see review bv lZurbuche'nl[20Q7h . 

The most remarkable exceptions are the numerical 
studies made in small cartesian domains which in- 
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elude the photospheric layers in the domain, the typ- 
ical numeric al domain size being of order of a few 
tens of Mm (iRempel et all [20091: iLeenaarts et al.l 120091: 



iMartfnez-Svkora et all 120081: 'Vogler et al .l 120051 among 
others; see also review bylwide mever- B6hm_eLalJ|2009). 
These studies are unable, though, to capture the slowly 
varying magnetic field and solar wind's properties at a 
global scale. 

The structure of the magnetic field inside the sun is 
unreachable by direct observation and/or measurement. 
Helioseismology techniques have been used, though, to 
deduce some of the properties of the ma gnetic field in 
the sun fe osovichev 2006|; IKosovichev fc Duvalll 120061 : 
lAntia et al.i ,2000, ,2003) . On the other hand, dynamo 
models are able to reproduce photospheric observables 
in greater detail. 

Currently our understanding of the inner solar mag- 
netism, its structure, evolution and origin rely on multi- 
D models of solar dynamo. This "fluid" dynamo is at the 
origin of the intense magnetic activity of the sun. In par- 
ticular, the solar interface dyn amo paradigm has received 
much attention ()ParkeHll993f ). This model assumes that 
the locations of generation of the toroidal and poloidal 
global field are separated, with the tachocline playing a 
central role in organising toroidal field. 

The classical explanation for the cyclic activity of the 
large scale magnetic field is that a dynamo process acts 
in the solar interior to regenerate the three components 
of the magnetic field and sustain them against ohmic 
dissipation. The inductive action of the complex fluid 
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motions would thus be responsible for the vigorous re- 
generation of magnetic fields a nd for its nonhnear evo- 
lution in the solar interior (see iCharbonneaul 120101 and 
iMieschI 120051 for recent reviews on the subject). Un- 
derstanding how these complex physical processes op- 
erating in the solar turbulent plasma non-linearly in- 
teract is very challenging. One successful and power- 
ful approach is to rely on multi-dimensional magnetohy- 
drodynamic (MHD) simulations. In this context, two 
types of numerical experiments have been performed 
since the 70's: kinematic mean-field axisymmetric dy- 
namo niodBls_wliich_^olvBOi^^ 

tion (ISteenbeck fc Krausd fToel IRobertslflQTl IStbaUOTGl: 
iKrause fc Radleiill980D and full 3D global mode ls which 
explicitly solve the full set of MHD equ ations (iGilmanI 
[1981 lGlatzmaieill985t ICatta nco 1999; Br un et al.ll200i r 

Clearly, both approaches are complementary and are 
needed to better understand the magnetic solar activ- 
ity. Recent progress have been made with 3D numeri- 
cal models of magnetic stars. Large-scale magnetic cy- 
cles are starting to be foun d in simulatioii s of th e Sun, 
such as those performed bv Ghizaru et al.l (|2010( ). or of 
solar-like stars (Brown et al. 2011). However, butterfly 
diagrams similar to the solar observations — with the 
poloidal field field reversal happening when the toroidal 
field is at its maximum — are still difficult to repro- 
duce. Kinematic mean-field models and their associated 
simplifying assumptions have thus been used extensively 
to reproduce several features of the large-scale solar cy- 
cle. In particular the use of a differential rotation profile 
inferred from helioseismology associated with an alpha- 
effect (due to the helical turbulence of the stellar convec- 
tive env elope) antisymmetric with respect to the equator 
enabled i Charbonneau fc MacGregorl (|1997D to produce 
a solar-like butterfly diagram with the ingredients of the 
modern interface dynamo. A similar model will be con- 
sidered in this work to catch the large-scale behaviour of 
the inner solar magnetic field. 

Above the photosphere, a very complex structure 
comprising open magnetic flux-tubes (coronal holes) 
and magnetic loops with different length-scales arises 
(magnetic carpet canopy, coronal loops) (jStixl 120021 : 
lAschwandenl l200"5l and references therein). The wealth 
of ground and space based observational data gathered 
in the latest years has been bringing up a great in- 
sight on the magnetic structure in the lower atmospheric 
layers (chromosphere and lower corona). The proper- 
ties of the magnetic field in the (optically thin) higher 
corona are harder to determine observationally, though. 
Space-borne in-situ measurements may introduce further 
constraints about how the solar wind properties con- 
nect to the coronal structures. Remarkable examples of 
these were given by the Ulysses consec utive polar orbits 
(jMcComas et al.l[200l Ifssautier et al.l[200l [2008 ). The 
SoHO spacecraft further provided complementary in-situ 
and solar surface (and low coronal) observations. Future 
missions such as the Solar Orbiter should provide more 
refined data to complete the scenario. These should com- 
bining high resolution imagery and spectroscopic data, 
magnetogram and in-situ measurements of the wind's 
properties in an orbit with very close helio-synchronous 
passes and varying orbital inclinations. 

The solar wind's outflow properties depend on the par- 
ticular geometry of the coronal field at each moment 



of the solar cycle. The amplitude and latitudinal dis- 
tribution of the solar wind velocity and mass flux de- 
pend on parameters such as the positions of the wind 
sources at the surface, the local magnetic field strength 
and the expansion factor of each particular coronal flux- 
tube the wind flows along. The details about the coro- 
nal heating mechanisms, most notably the amplitude 
and loc ation of energ y dissi pation, are also of impor- 
tance. [Leer fc Holzeii ([1980f ) pointed out that heating 
below or above the critical sonic point may produce dif- 
ferent effects on the wind velocity, and that high speed 
winds require ene r gy de position in the supersonic region. 
iHansteen fc Leeil ( 1995[ ) showed that the mass flux does 
not depend strongly on the mode and location of en- 
ergy deposition, but rather on the amplitude of the en- 
ergy flux. Note also that the heat flux profile could be 
related to the i nagne tic field's am plitud e and geometry 
(iCranmer et al.l 120071 : [Pinto et al.l 120091 ). The resulting 
angular momentum losses and wind's braking torque (ap- 
plied on the sun) depend on wind's velocity and mass 
flux together with the (time varying) geometry of the 
Alfven surface (the surface at which the wind's velocity 
equals the Alfven velocity; see deflnitions in f}5]). Some 
authors have studied the influence of the topology of 
coronal field in the properties of stellar winds and re- 
sulting braking torques, but mostly by using simplified 
config urations (e.g. dipolar versus quadrupolar field, as 
in Ma tt fc Pud ritz 2008). Others have furthermore con- 
sidered how the presence of strong magnetic spots (lo- 
calised magnetic fiux enhancements) specifically affects 
the angular mome ntum loss rate fe.g. lAibeo et al.ir2007l : 
ICohen et al. [ [2009l) . 

More detailed studies of the global solar magnetic 
structure usually adopt surface magnetogram data for 
the radial component of the field as a lower bound- 
ary. The atmospheric magnetic field's geometry is de- 
duced using potential field extrapolation techniques. 
The coronal and heliospheric hydromagnetic condi- 
tions are then ded uced eith er by using a set of scmi- 
empi r ical relations (|Arge fc Pizzo 2000; Luhmann et aO 
[200l IWang fc S hcclcv 1992: Schriivcr fc DcRosa 2003), 
or b y finding stable MHD solutions for the mapped fields 
(e.g. lUsmanov[[l99l[Huet al.ll2008D . The usage of mag- 
netogram data as lower boundary conditions encounters 
a few difficulties. Namely, only the line-of-sight com- 
ponent of the field is available, making the estimation 
of the radial components difficult at the polar regions. 
The standard PFSS technique (standing for "potential 
field source surface") further assumes the surface field 
is strictly radial. One still needs to assume an outer 
boundary condition which imposes the magnetic field to 
be purely radial as well at a source- surface (typically, a 
spherical surface placed at about r = 2.5 R©). This con- 
dition emulates the field-line opening usually caused by 
the solar wind fiow, while keeping the global field at its 
lowest (current-free) energy state. These studies essen- 
tially provide quasi-stationary predictions of the coronal 
fields at a given moment. Despite all the simplifications 
and assumptions this family of models has been success- 
ful in predicting qualitatively the coronal magnetic field's 
topology (observed, for example, in white light during 
solar eclipses) and became a standard to which other 
types of model can be compared. Their MHD counter- 
parts are limited in different ways. Beyond the much 
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increased computational power they require, they rely 
on phenomenological assumptions for the heat transport 
in the corona (as the heating processes are still under 
debate). These models show, nevertheless, some quan- 
titative differences when compared to the PFSS ones. 
The relative sizes of the streamers can be different (usu- 
ally more elongated) and the heliospheric magnetic field 
tends to be more unif orm (as it is known to be). See, for 
example, iRilev et al.l ([2006) for a comparison between 
these two families of models. A third type of models — 
nonlinear force-free field (NLFFF) models — places it- 
self somehow in between the two previously described, 
in the sense that it tries to add some complexity to the 
PFSS scenario while avoiding the issues found by global 
MHD computations. These require, though, good qual- 
ity vector magnetogram data at the surface of the sun 
and still find some difficulty in accurately predicting real 
solar features (DcRosa et aL 2009). 

In this paper, we investigate the influence of the cyclic 
evolution of the large-scale magnetic field produced by 
the solar dynamo on the solar wind properties. We chose 
not to perform this study based on surface magnetogram 
magnetic field data (as in Wang et al. 2003; Hu ct al. 
I2008t iLuhmann et al.ll2009l for example), but rather pur- 
sue a diff'erent path: our source fields are the result 
of well-tested and validated kinematic dynamo models 
(jJouve et al.ll2008f) . An isothermal MHD model of the 
solar wind and corona is used to produce a temporal se- 
quence of steady-states spanning a complete activity cy- 
cle. The focus here is on estimating the wind's velocity, 
mass and angular momentum fiux spatial profiles as they 
vary during an activity cycle in response to the variations 
of the magnetic field's topology, rather than reproducing 
a particular solar cycle. The variability of global proper- 
ties such as the sun's mass loss rate and wind's breaking 
torque are studied in regards to the activity cycle. Also, 
we expect to gain a deeper insight on the connections 
between the sub-surface and coronal physical processes 
by proceeding this way. 

The remainder of this manuscript is organised as fol- 
lows: the methods and numerical codes used are pre- 
sented in Sj2j The results are described thoroughly in [J3] 
(coronal magnetic field evolution) , gl (solar wind speed) 
and fJS] (mass and angular momentum flux, breaking 
torque) . A discussion follows in ^ and a brief summary 
in §3 

2. COUPLING SOLAR DYNAMO AND WIND MODELS 

We used two 2. 5D axisymmetric MH D codes. The flrst 
one — STELEM pouve fc Brun|[2007l and gUI]) — com- 
putes mean-fleld kinematic MHD dynamo solutions given 
the meridional circulation and differential rotation pro- 
flies in the solar convection zone. A Babcock-Leighton 
or an a source terms generat e the poloidal ma gnetic 
fleld. The second code — DIP ()Grappin et al.ll2000 . and 
i j2.2p — computes the temporal evolution of an MHD 
solar corona with a self-consistent wind. In the latter, 
the magnetic fleld separates into two components: an 
imposed stationary external potential fleld B*^ (whose 
sources lie within the sun), and an induced component 
produced by the flows within the numerical domain. The 
dynamo magnetic fields produced by STELEM match a 
potential field at the surface of the sun. This potential 
source field can then be directly transmitted to DIP as 



the external component of the coronal magnetic field. 
Further details about both numerical codes are given 
hereafter, in !j2.1l and Details about the coupling of 
the two codes are given in §2.31 




Figure 1. Representation of the numerical domains of STELEM 
(left) and DIP (right) at the same instant. The colour scale in 
the left figure shows the toroidal magnetic field B^f, in the convec- 
tion zone. Black lines are (poloidal) magnetic field lines (dotted 
and continuous lines indicates CW and CCW field-lines, respec- 
tively). The white lines in the right figure are magnetic field lines, 
while the colour-scale represents both the wind's velocity in units 
of Mach number and the open field's polarity (blue is negative, red 
is positive). 



2.1. STELEM: a representative dynamo solution 

To investigate the global solar cycle features produced 
by dynamo models, we start from the hydromagnetic in- 
duction equation, governing the evolution of the mag- 
netic field B in response to advection by a flow field U 
and resistive dissipation. 

r)B 

_ = Vx(UxB)-Vx (r;V x B) 

As we are working in the framework of mean-field the- 
ory, we express both magnetic and velocity fields as a sum 
of large-scale (that will correspond to the mean field) and 
small-scale (associated with fiuid turbulence) contribu- 
tions. Averaging over some suitably chosen intermediate 
scale makes it possible to write two distinct induction 
equations for the mean and the fiuctuating parts of the 
magnetic field. 

A closure relation is then used to express the mean elec- 
tromotive force in terms of the mean magnetic field, lead- 
ing to a simplified mean- field equation (Moffatt. 1978) . 
The prescribed mean velocity field will here only con- 
sist in its longitudinal component (the solar differential 
rotation) and the magnetic diffusivity is assumed to be 
constant. The source term for the poloidal field is linked 
to the turbulent helical motions within the convection 
zone. We thus obtain a simple af2 dynamo model with 
constant diffusivity. 

Working in spherical coordinates and under the as- 
sumption of axisymmetry, we write the total mean mag- 
netic field B as 

B(r, 0, t) = V X (A^(r, 9, t)e^) + B^{r, 9, t)e^ 
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Reintroducing this poloidal/ toroidal decomposition of 
the field in the mean-field induction equation, we get 
two coupled partial differential equations, one involving 
the poloidal potential and the other concerning the 
toroidal field B, 



value of fleq = 460 nHz, in agreement with observations). 
The spatial resolution is 128 x 128. 







dt 
dt 



Cc^aB^ + (V^ 



(1) 



(2) 



where w = rainO, f2 is the differential rotation and a is 
the Qf-effect. 

All these quantities are now dimensionless, thanks 
to the presence of the two Reynolds numbers Co = 
^eqR%/rit, and Ca — ctoRQ/rjt, rjt being the turbulent 
magnetic diffusivity, ileg a measure of the rotation rate 
at the equator and ao a measure of the intensity of the 
alpha-effect. The product of these two numbers gives us 
the dynamo number which measures the efficiency of the 
source terms to make the magnetic energy grow in time 
against Ohmic diffusion. 

Equations ([T]) and ^ are solved in an annular merid- 
ional cut with the colatitude 6 g [0, tt] and the radius 
r S [0.65, 1] R© i.e from slightly below the tachocline 
(e.g. r = 0.7 Rq) up to the solar surface (see Figure 
[11 left panel). At = and 9 = tt boundaries, both 
A^ and B^ are set to 0. At r = 0.65 Rq, we compute 
a perfect conductor condition. At the upper boundary, 
we smoothly match our solution to an external potential 
field, i.e. we have vacuum for r > Rq. 

The STELEM (STellar ELEMents) code uses a finite 
eleme nt method in space and a third o rder scheme in 
time (|Burnettlll987llJouve fc Brunl[200l . 

The principle of the finite element method is to look for 
solutions of the weak formulation of the equations, here 
equations ([T]) and ([2|). Those solutions are taken to be 
linear combinations of well-chosen trial functions. In our 
case, they are Lagrange polynomials of degree 1 and the 
elements are rectangles in the (— cos 9, r) plane. Apply- 
ing this spatial method to our equations results in a sys- 
tem of first order ODEs in time governing the evolution 
of the coefficients of the linear combinatio ns. The tem- 
oral scheme that we use is adapted from iS palart et al.l 
1991h . It is similar to a Runge-Kutta 3 method and is 
thus explicit and of order 3. 

The STELEM code is here used to produce a cyclic 
dynamo field within the model convection zone, whose 
values at the top of our domain will be reintroduced in 
the DIP code described below. We chose a simple afl 
dynamo model to produce our representative solution, 
with physical source terms consistent with observations. 
The rotation profile is deduced from helioseismic inver- 
sions and the a-effect is antisymmetric with respect to 
the equator and positive in the Northern hemisphere (in 
agreement with the preferred handedness of the turbu- 
lent helical convective motions) . This is similar to case 
A of pouve et al.ll2008( ). The following values were used 
for the various parameters: Ca = 0.385 (critical value 
for which dynamo action just starts to compete against 
Ohmic diffusion), Cn = 1.4 x 10^ (corresponding to a 




Time 

Figure 2. Time-latitude diagram of the surface field. The upper 
panel shows Bg while the lower panel shows Br- Red (blue) colours 
indicate positive (negative) values of the field. The white dashed 
lines span the time interval which was chosen to be incorporated 
in the DIP code to compute the solar wind evolution and mark the 
instants shown in the following figures. 

Figure [2] shows the temporal evolution of both the lati- 
tudinal and the radial field at the surface, at all latitudes. 
We have here the evidence of a cyclic magnetic field: the 
surface poloidal field changes sign regularly and at dif- 
ferent times depending on the latitude of interest. This 
simple model does not aim at reproducing the solar dy- 
namo in its finest details. The goal here is rather to 
assess the infiuence of a cyclic poloidal background field 
(having the same general properties as that in the sun) 
on the wind's properties. The dashed white lines labelled 
ti and ig in the figure indicate the beginning and the end 
of the particular time interval used to compute the so- 
lar wind evolution. This time interval spans an entire 
cycle period and is roughly located between two activity 



2.2. DIP: wind model 

The DIP code is a 2.5 D axisymmetric model of the 
solar corona obeying the compressible MHD equations 
for a one-fluid, isothermal and fully ionised plasma. The 
continuity and momentum equations are 

dtp + V ■ pu = (3) 

5tU-H(u- V)u = 

VP J X B , 
\ g + j/V^u . 

We set 7 = 1 and a uniform temperature T = 1.3 MK 
in all the corona. The gas pressure is deduced from the 
equation of state 

P = —pksT , 
niH 

valid for a fully ionised hydrogen plasma. The isother- 
mal approximation discards a complete treatment of 
the energy fluxes in the corona. But we should note 
firstly that the corona is nearly isothermal (at least 
in the first 12 Rq). Secondly, assuming 7 = 1 is a 
proxy to the actual thermal state of the coronal plasma. 
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summing-up the combined effects of the thermal con- 
duction and the (still debated) heating source. The 
choice of a particular value for the coronal tempera- 
ture is somewhat arbitrary, though. The temperature 
T = 1.3 MK chosen here is justified both by represent- 
ing an averag e coronal value in the first ~ 10 R0 (see , 
e.g, results bvlHansteen fc Leeiiri 995: End eve et al.il2003l : 
iGuhathakurta et al.ll2006t iPinto et al...2009r among oth- 
ers) and empirically as it produces correct wind solu- 
tions. As discussed later in the text, the wind's mass flux 
depends both on the coronal temperature and on wind 
geometrical expansion factors. We focus here on the geo- 
metrical effects linked to the activity cycle and keep our 
reference temperature fixed throughout the simulations. 

The magnetic field B decomposes into a potential ex- 
ternal component B" and on a component b induced 
by the flows. The former only has poloidal compo- 
nents (r, 9), while the latter comprises both poloidal and 
toroidal components. The total magnetic field is then 
defined as 

B = B" + hpoi + . (4) 

Furthermore, hpoi = V x We integrate the evolution 
equation for the magnetic potential # 

9t* = UrBe - Brug + r/V^* . (5) 

The poloidal components of b are then computed as 

r 

$ 

hg = . 

r 

The azimuthal component of the induced field comes di- 
rectly from the induction equation (in the direction 
only), that is 

dth^^ [Vx (uxB)+ryV2B]_^ . 

This method guarantees that the solenoidal condition V • 
B = is satisfied at all times. This papers focus only on 
the poloidal components of the coronal field, though (as 
explained hereafter, in i j2.3p . 

The diffusive terms are adapted so that grid scale (AZ) 
fluctuations are correctly damped. The kinematic vis- 
cosity is defined 'as v = Vfj (AZ/AZq)^, typically with 

z/Q = 2 X 10" cm2-s-i and 0.01 < {Al/Alof < 10. 
The simulations presented in this paper were performed 
with a vanishing 77, so as to approach as much as pos- 
sible the limit of ideal non- resistive MHD. Additional 
diffusive terms are used. These are implicit numerical 
filters (applied over u, b and $) which dissipate mostly 
at the grid scale a nd minimi se the dissipation of large 
scales fluctuations (jLeldllQQ^ . This flltering scheme al- 
lows the diffusive parameters u and j] to be lowered while 
avoiding spurious G ibbs fluctuations, and all ow for lower 
mid-scale damping (jGrappin fc Leora^l2001|) . Note that 
actual kinetic dissipation should happen at scales much 
smaller than the grid size, anyway. 

We used a 512^ grid which is uniform in latitude and 
non-uniform in radius. The grid's cells radial extent is 
Sr — 6.5 X 10"'^ Rq at the lower boundary and Sr = 
1.0 X 10~^ R© at the upper boundary. 

Both the upper and the lower numerical boundaries 
(respectively at r = 15 and 1.01 Rq) are transparent. 



i.e, mass can flow through them and waves are not 
spuriously reflected there (see Figure [U right panel). 
This is achieved by writing the MHD equations in 
their characteristic form. In simple terms, this con- 
sists in projecting the system of equations in terms of 
the primitive quantities {p, U,B} into an equivalent 
system defined in terms of characteristic variables Lj 
(iThoinpsonlllQSTHVanaiakshi et al.lll989l:lPoinsot fc Leld 
II992I : iRoe fc Balsaralll996[ ). The evolution equation for 
these variables explicitly define the time-evolution of the 
MHD system in terms of wave-modes propagating up- 
wards and downwards, or equivalently, incoming and out- 
going modes. We must then constraint all the incoming 
modes, and let the others free. Plasma is free to flow 
through and its properties to vary in time, the actual 
boundary values at each moment depending on the cur- 
rent state of the system. The mass fluxes and velocities 
(or their gradients) are, therefore, not arbitrarily set at 
the numerical boundaries. This is a critical feature: the 
actual solar wind velocities and mass flux at each point 
of the surface will automatically match each unique tran- 
sonic and transalfvenic wind solution. There is no need 
to set in advance arbitrary values for the density, ve- 
locity or mass flux at the numerical boundaries. This 
means no spurious boundary layers will form there and 
that the resulting mass fluxes will not be artiflcially con- 
strained. We note that the resulting surface density and 
mass fiux do not vary by large factors (not more than 
^ 10% in all latitudinal domain, and less than 6% within 
open field regions at all moments of the cycle). Further- 
more, this type of boundary conditions better describes 
dynamic finite-frequency phenomena th an any type of 
non-t ransparent (i.e, "rigid") conditions (jGrappin et al.l 
l2008h . 

The solar wind develops into a stable transonic and 
transalfenic solution in the open field regions after a pres- 
sure (sonic) perturbation is applied at the outer bound- 
ary. This perturbation propagates inwards until both 
the sonic and alfvenic surfaces appear (and are enclosed 
in) within the domain. From then on, the solar wind 
profiles evolve towar ds a unique stable state (Velli 199^ 
iGrappin et al.l Il997| ) . Note that the frontiers between 
open (coronal holes) and closed flux regions (streamers) 
are not set in advance, but they are rather a result of 
the competition between the magnetic tension and the 
dynamical pressure in the wind flow. In other words, the 
solar wind settles where the magnetic field is incapable 
of containing the outward flow and forces the magnetic 
field to align with the radial direction. 

Other aspects of the numerical model are more thor- 
oughly discussed in IGrappin et al.l (|2000D . 

2.3. Coupling method 

We used a t ime-series (t) ,B a (t)} derived from a 
STELEM run ()Jouve fc Brunll2007D which correctly pro- 
duces a cyclic behaviour of the poloidal magnetic field. 
We used a simple afl dynamo model with source terms 
for both the toroidal field and the potential A^. 

The time-series was scaled in order to span an 11-year 
period, sampled with a time step of 6 months. The verti- 
cal dashed white lines labelled ti and te in Figure [5] show 
the beginning and end of the time-series (^i — and 
^6 = 11 yr). These are the moments when the poloidal 
field is at its simplest configuration (i.e, lowest multipo- 
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lar order). For simplicity, we will call these monients 
"activity minima" , while the magnetic field's polarity in- 
version phase (between lines t2 and is in the same figure) 
will be called "activity maximum" hereafter. The time- 
series comprises 22 equally spaced samples, but we will 
refer to a subset of 6 samples {ti to t^) for illustrative 
purposes throughout the text. We let the solar wind f ully 
develop in the domain and reach a steady state (see ^2.2\ 
for further details on the formation of the wind.) The 
amplitude of the external magnetic field B*^ was scaled 
so that the total field would match coronal amplitudes 
(the scaling factor being kept constant during the whole 
cycle). Then, for each consecutive sample, we substi- 
tuted directly i — ^ B^^^. The system was again let 
free to relax and attain a new steady state for each itera- 
tion i. The stability of the relaxation method was tested 
so that we could be sure that each relaxed steady-state 
solution did not depend on the history. That is, we found 
no hysteresis when cycling back and forth through dif- 
ferent states (given that we let the system relax at each 
stage) . Different permutations of the original time-series 
i = 0,1, . . . , N also lead to the exact same steady-states 
as long as each solution's set of parameters was kept the 
same. The wind solution for each instant of time i de- 
pends only on the corresponding B^. 

Some of the runs were performed at different grid res- 
olutions, and numerical convergence was verified. 

Both numerical codes are time-dependent in nature, 
but the result of the whole coupling procedure is a se- 
quence of steady-state solutions. Our procedure gener- 
ates a map of the poloidal coronal magnetic field and 
wind fiow during an activity cycle rather than the dy- 
namical evolution of a particular event. 

Toroidal fields and flows will not be considered 
in this paper. Concerning rotation, the sun is a 
slow magnetic rotator, and therefore the magneto- 
rotational effects on the polo idal wind flow are negligible 
(|Belcher fc MacGregod[l976 h. 

3. CORONAL MAGNETIC FIELD 

The temporal evolution of the corona and solar wind in 
response to the dynamo field variations is shown in Fig- 
ure [21 Only the first 4 Rq of the northern hemisphere 
are displayed at six different instants of the activity cy- 
cle, corresponding to the ti, t2, ■ ■ ■ , lines in Figure [2l 
That is: t = 0, 3.3, 3.6, 3.8, 4.4, 11 years. The colour- 
scale represents the solar wind's velocity projected on the 
unit magnetic field vector in units of Mach number, that 
is V • B/(cs||B||). Orange/yellow and blue/green shades 
therefore trace different B-field polarities in the open 
field regions (respectively, u • B positive and negative). 
The sharp transitions between positive and negative po- 
larities in this figure outline current sheets (note that the 
wind fiow does not change sign across these transitions, 
but Br does). 

Some elements and characteristics are observed consis- 
tently throughout the whole cycle. Higher concentrations 
of magnetic flux at the surface (or equivalently, of cur- 
rent below the surface) translate into coronal loop arcade 
systems. Strong flux concentrations appearing in coronal 
holes shape up as helmet streamers, and end in a cur- 
rent sheet which extends outwards. Smaller flux concen- 
trations embedded in unipolar flux regions form nearly 
symmetric bipolar structures with no current sheet (as 



the "giant polar plum es" in iPinto et al.l 120101 and the 
"pseudo-streamers" in iWang et al.ll2007n . Loop struc- 
tures placed inside stagnant zones (i.e with no wind fiow, 
as inside large streamers) mostly maintain their potential 
field configuration. 

We start the computation at the moment when the 
magnetic structure of the corona is at its simplest (Fig- 
ure [3l first panel, ti ; also first white line in Figure [2]) . 
One and only large equatorial streamer extends from the 
surface up to nearly 5 Rq, where the heliospheric cur- 
rent sheet starts. The frontiers between the streamer 
and the coronal holes cross the surface at latitudes —60° 
and -|-60°. The letter A indicates the latitudinal extent 
of the streamer. The open magnetic field has positive 
polarity in the northern coronal hole and negative polar- 
ity in the southern coronal hole. The streamer itself is 
divided into four magnetic connectivity regions around 
one X-type null point. The null point itself is located 
over the equator at r — 1.4 Rq. The four connectiv- 
ity regions mentioned above are then the group of small 
equatorial loop arcades place below the null point, the 
group of larger equatorial arcades above the null and 
filling up most of the streamer, and the two groups of ar- 
cades to the north and south of the null (note that Figure 
[3] only shows one hemisphere, and that the system is sym- 
metric with respect to the equator). As the solar cycle 
starts moving away from the minimum, new fiux con- 
centrations emerge at mid-latitudes (showing up as new 
groups of coronal loops inside the equatorial streamer). 
These new structures slowly migrate polewards, attain- 
ing the streamer boundaries at about t = 2.5 yr. The 
equatorial streamer is disrupted at this point. A fraction 
of the magnetic fiux will remain in the equatorial region, 
forming a smaller streamer (Figure [3l second panel, ^2)- 
The rest of the magnetic fiux reconnects and forms new 
pl ume/ pseudo- stream er structures at mid- latitudes (as 
in IPinto et al.l l20Tol and iWang et al.l [20071 ) . The letter 
B indicates the pseudo-streamer position. New coronal 
holes now appear at low latitudes. The polewards pro- 
gression continues at a steady pace opening up its way 
by reconnecting with the open magnetic fiux. At about 
t = 3.5 yr one of the magnetic arcades of each of the 
newly formed pseudo-streamer breaks and quickly opens 
up (Figure [21 third panel, ^3, letter C). As a result, a 
new coronal hole with inverse magnetic polarity rapidly 
forms. This coronal hole will grow wider and fill up all 
the polar regions as the polarity inversions proceeds. The 
previous coronal hole's field then closes down (Figure [H 
fourth panel, ^4), and ends up disappearing below the 
surface (Figure [31 fifth panel, t^). The corresponding 
coronal arcades (the ones closing down near the poles) 
are indicated by the letter D in the figure. At t = 11 yr 
all traces of the previous coronal hole magnetic field have 
disappeared (Figure[3l last panel, t^). The system is back 
to a state very close to its original state, but with the po- 
larity of the magnetic field reversed (the dynamo model 
used here produces very regular and symmetric cycles). 

Note that the polarity reversal happens quickly in the 
corona, even if the underlying B^'-field evolves slowly. 
The wind flow is responsible for the quick "opening-up" 
of field lines as the magnetic flux concentrations evolve 
slowly at the surface, and ultimately for the change of 
connectivity between contiguous regions. Furthermore, 
at the polar axis, the magnetic field's inversion occurs at 
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Figure 3. Snapshots of the evolution of the corona during the solar cycle (only the first ~ 4 Rq and northern hemisphere are shown) at 
the instants ti to shown in Figure [2] White lines are magnetic field lines. The colorscale represents the quantity ^rjg^i that is, the 

wind flow velocity projected onto the signed magnetic field in units of Mach number. This quantity traces the B-field's polarity in the open 
field regions. Red/orange means positive polarity, while green/blue means negative polarity. The large arrowheads show the local B-field 
orientation. The grey contour shows the sonic surface. The letters A, B, C and D indicate the positions of particular magnetic structures 
which we refer to in the text. 
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Figure 4. Top panel: Polar magnetic field as a function of radius 
at different moments of the cycle. The instants represented are 
the same as in the panels in Figure |3](ti to ig) plus an additional 
t^h = 5.5 yr for completeness. Bottom panel: B^r^ at the north 
pole as a function of time. B is in units of G and r is in units of 
R0. Each curve corresponds to a different height. The factor 
accounts for the field's decay due to a purely radial expansion (note 
that Brr'^ decays faster in the lower part of the domain, but not 
above) . The polarity inversion at the surface is delayed with respect 
to higher coronal heights. This is due to the slowly progressing 
closing down of polar fields, as seen in Figure [3] 

different times at different fieigfits. To better describe 
tfiis property, the top panel in Figure 2] sliows radial cuts 
of the polar magnetic field at the same moments as those 
shown in Figure[3l The bottom panel shows the temporal 
evolution of the polar field at different heights. The Br 
sign-switch happens first at higher altitudes and proceeds 
downwards. The downwards progression of the reversal 
of Br is quick {6t of order of a few days) between r = 
15 R© and r = 2 Kq, but from there on it slows down 
considerably. The delay is « 6 months between r — 
2 R0 and r = 1.6 R©, and about 1 yr between r = 
1.6 R© and r = 1.3 R©. The overall delay (between 
top and bottom) is about 4 yr. This delay in the lower 
coronal layers is due to the slow disappearance of the 
coronal arcades near the poles, as shown in the last three 
panels of Figure |3l 

Figure [5] shows Br as a function of Q at different al- 
titudes and at different moments of the cycle. Faraway 
from the surface, B « Bre,r and is nearly independent 
of the latitude, except in the vicinity of a current sheet. 
Close to the surface (beneath r « 2 — 3 R©), the mag- 
netic field organises in much more complex ways. The 
potential component of the magnetic field dominates in 
the lower layers of the corona whereas the induced com- 
ponent largely dominates above (as the wind flow be- 
comes stronger and approaches its asymptotic velocity). 
From there on, the solar wind flow and the total magnetic 
field essentially align with the radial direction. 
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Figure 5. Brr^ as a function of latitude at different heights, and 
at different moments of the cycle. The factor accounts for the 
field's decay due to a purely radial expansion, as in the previous 
figure. B is in units of G and r is in units of Rq . The instants 
represented are the same as in the panels in Figure [3] (ti to tg) 
plus an additional t^i, = 5.5 yr for completeness. Faraway from the 
sun, the radial magnetic field is mostly uniform in latitude (except 
around current sheets, where it changes sign), independently of the 
complexity of the surface field. 

The multiple current sheets shown in the intermediate 
panels in Figure El may be interpreted as a highly warped 
current sheet in the real non-axisymmetric corona. 

At the activity minimum the open fiux is restricted 
to large polar coronal holes (about 30° in latitude in 
each hemisphere). This open fiux will eventually fill all 
the available space at greater heights (spanning 180° in 
latitude, from pole to pole), as the streamer thins out; 
above r « 3 R© all the magnetic field is open. Closer to 
the maximum, the open fiux sources are more spread in 
latitude but cover a smaller latitudinal extent altogether 
(about 5° in latitude in each hemisphere). These multiple 
thin coronal holes will nevertheless grow with height, and 
will end up filling up all available space. Above r w 
2.5 R© all the magnetic field is open. In other words, the 
average flux-tube expansion factors will be much higher 
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at the maximum than at the minimum. 

4. SOLAR WIND SPEED 
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Figure 6. Solar wind speed at the outer boundary during ttie cy- 
cle as a function of time and latitude. Most of the time, faster wind 
outflows occupy a large latitudinal extent ranging from the poles 
to latitudes as low as ±20° . The exception is the polarity reversal 
phase (t £ [2,5] yr). Two low-latitude wind channels appear at 
about t = 3_yr, corresponding to the newly formed coronal holes 
(see Figure [3] at (3). 

We focus now on the variation of the solar wind veloc- 
ity during the solar cycle, and on how it is distributed in 
latitude. The evolution of the coronal magnetic topology 
(f}3|) has a direct influence on the size and distribution of 
wind sources at the surface of sun, via changes in posi- 
tion and width. The local magnetic pressure dominates 
over the wind's dynamical pressure at the lowest layers, 
and the local B-field's amplitude and inclination mostly 
determine whether a given field line is open or closed. 
Conversely, the latitudinal distribution of the solar wind 
faraway from the sun cannot be trivially predicted from 
the magnetic field's configuration at the surface. The so- 
lar wind flows accelerates along open field lines. A fiuid 
element of cross-section Aq at the surface will accelerate 
along a magnetic flux-tube with cross-section A (r) . The 
final velocity profile depends on the flux-tube's expan- 
sion factor . The expansion factor itself results from 
the competition between the wind's dynamical pressure 
and the magnetic pressure in the corona. 

The terminal wind speed and flux tube expansion 
factor are inversely correlated at all latitudes and at 
all times, agreeing with the well es tablished Wang- 
Sheeley-Arge semi-empiri cal relation (jWang fc SheelevI 
119901: fArge fc Piz^ I2OOO0 . In the initial acceleration 
phase (below r ~ 3 R©) the expansion factors vary con- 
siderably, both in latitude and radius. At larger radii, 
though, the multiple and initially independent wind flows 
merge into a bulk spherical outflow {A(r)/Ao oc r~^, 
independently of the latitude). Latitudinal inhomo- 
geneities still subsist in this flow, being a result of the 
evolution (acceleration) of each parcel of wind flow along 
its path starting from the sun's surface. 



Figure [S] shows the latitudinal distribution of the solar 
wind speed at r = 15 Rq (the domain's outer bound- 
ary) during the solar cycle. The wind velocity vs. lati- 
tude diagram shows g ood qualitative agreeme nt with the 
predictions made by IWang fc SheelevI (j2006[) matching 
ULYSSES in situ measurements and the more recent 
multi-station I PS (interplanetary scin tillation) coronal 
observations bv lTokumaru et aTl ()201Cl[ ). The separation 
between fast and slow wind is well visible in the figure. 
At the activity minima, fast solar wind originates essen- 
tially from high latitude regions, while the slow wind 
fiows mostly closer to the streamer frontiers, at lower 
latitudes. As the activity cycle progresses from the min- 
imum to the maximum, the slow wind expands over to- 
wards the poles and takes over most of the latitudinal 
domain. On the declining phase of the activity cycle, the 
fast wind recovers the polar regions and progressively ex- 
tends towards lower latitudes, restraining the slow wind 
flow to the equatorial region. Some irregularities appear 
above this otherwise too simple scenario. Most remark- 
ably, two short-lived wind channels appear at low lati- 
tudes between i = 3 yr and t = A yr, that is, during the 
polarity inversion. This can be seen in Figure |3] (pan- 
els t2, t2, ^3); note how the blue shaded wind channel 
appears and evolves. These flows originate in the newly 
formed mid-latitude coronal holes. They follow the corre- 
sponding magnetic flux-tubes, bending over from 50° 
at the surface down to ~ 30° at the outer boundary. 
Although they seem to fade away quickly in Figure IHl 
they actually last till the declining phase of the activity 
cycle, but the wind flowing within these coronal holes 
slows down. This slowing down is due to the increas- 
ingly higher expansion factor for the coronal hole, as can 
be seen in Figure [31 Note how the blue shaded coronal 
hole expands to fill the whole hemisphere faraway from 
the sun while its latitudinal extent at the surface remains 
approximately constant (from t^ to ^5). This continues 
while the polar closed-field regions progressively disap- 
pears; all open flux will merge into a wide polar coronal 
hole afterwards (ig)- Also, the closing down of magnetic 
flux near the poles (during the cycle's decay phase, at 
about t — 7 yr) is related to the appearance of fast wind 
flows close to the polar axis. These correspond to newly 
formed thin polar flux-tubes expanding almost radially, 
which will also merge afterwards into the wide polar coro- 
nal hole. 

The wind's velocity values presented here are expected 
to be lower than the values measured in situ near the 
Earth's orbit. The reader should note that, on one hand, 
our numerical domain extends only up to 15 R©, and 
that the wind flow has not yet reached its asymptotic ve- 
locity at this height. Nevertheless, the relative variations 
of Vr in latitude should not change considerably. The 
wind undergoes a purely spherical expansion between 
15 R0 and 1 AU. On the other hand, this does not fully 
accounts for the low velocities values found, though. We 
could expect velocities in the order of 500—600 km/s from 
this model if the domain extended up to 1 AU. Faster 
wind velocities require a complete and consistent ener- 
getic treatment, which is bey o nd th e scope of the current 
paper. Besides, iPinto etall ()2009[ ) show that the anti- 
correlation between terminal wind velocity and flux-tube 
expansion factor is still verified in self-consistent non- 
isothermal cases. The only strong constraint over the 
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domain's radial extent regarding the physical correctness 
of the numerical model is that it has to completely con- 
tain all critical surfaces (sonic and alfvenic; cf. 
which it always does. 

5. MASS FLUX, MOMENTUM FLUX AND MAGNETIC 
BRAKING TORQUE 
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Figure 7. Mass flux pVrr^ sin 9 in the meridional plane at t = yr 
(left) and at 4 = 3.3 yr, about the polarity inversion (right), that 
is, respectively, instants ti and t2. The factor r^sinS is due to 
the spherical expansion of a surface element normal to the radial 
direction. Outflows which originate at lower latitude dominate the 
global mass loss rate. 
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Figure 8. Mass loss rate M during the solar cycle. 

Figure [7] shows the radial mass flux associated with 
the solar wind in the corona at two different instants 
of the solar cycle. In both cases, the net mass outflow 
in the polar coronal holes is higher near the streamer 
boundaries than closer to the poles. At the maximum 
of activity, thin coronal holes appear also at low lati- 
tudes (as described previously). Despite their small lati- 
tudinal extent at the surface, the associated mass flux is 
important when compared to that in the polar coronal 
holes. For an outflow with a given latitudinal extent S0, 
the actual surface area it crosses is equal to 27rr^ sinOSB. 
Therefore, low latitude coronal hole are more prone to 
produce higher mass outflow rates. Arguably, this is a 
consequence of the axi-symmetrical nature of our model. 
Nevertheless, the real three-dimensional sun shows some 
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Figure 9. Alfven surfaces (black contours) at the same instants as 
in fig. [7] (instants ti and t2). White lines are magnetic field lines, 
and the colorscale represents the wind's poloidal Mach number. 



degree of axi-symmetry in the distribution of the stream- 
ers and coronal holes. These low-latitude wind streams 
correspond to the two channels visible in Figure [51 at 
about t — 3 yr. 
Figure [8] shows the total mass loss rate 

M^2'n:Rl[ pVrSm9d9 (6) 
Jo 

evaluated at the outer boundary of the numerical do- 
main. The mass loss rate evolves in par with the activity 
cycle. That is, M is maximal at about t = 3 yr (during 
the activity maximum) and minimal during the activity 
minima. The amplitude of the mass loss rate varies by 
a factor of about 1.6 in time, from 4.2 x 10"^'' Mq/yt 
at the activity minimum and 6.9 x 10^^* Mq/jt at the 
activity maximum. This trend supports the idea that 
the lower latitude outflows — which appear mainly close 
to the activity maximum — contribute with higher net 
mass outflow rates, as discussed in the previous para- 
graph. Note that the increase in M cannot be due to the 
variations in the velocity of the wind, as it is actually 
lower at almost all latitudes during the activity maxima 
(see Figure [6]) and contributes to lowering the net mass 
flux. On the other hand, the variations of p at the lower 
numerical boundary are small throughout all the cycle. 
The amplitude of the density fluctuations at the surface 
within open-field regions is always below 6% of its aver- 
age value (amounting to at most 10% if both open and 
closed-field regions are considered altogether). 

The way the open flux maps from the solar surface up 
the outer domain must therefore be the main cause for 
the variations found in the mass loss rate. 

The solar wind outflow carries angular momentum 
away from the sun. The specific angular momentum flux, 
magnetic braking torque and spin-down time-scale are 
deduced from the Alfven surface's geometry and solar 
rotation rate fig- 

The Alfven surface is the geometric locus where the 
wind velocity equals the Alfven speed ca — B/y^Anp. 
The Alfven radius ta is the cylindrical distance from 
the rotation axis to this surface. Classical wind theory 
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Figure 10. Mean Alfven radius {ta) (top), specific angular mo- 
mentum flux n(r^) (middle) and magnetic braking time-scale 
^isd = Jq/J (bottom) during the solar cycle. Jq is the sun's 
angular momentum IIGilman et al.|[l989l) : J is the wind's angular 
momentum loss rate. 



(jWeber fc Davislll967l and many others) states that the 
angular momentum balance problem can be simplified 
as follows. The plasma inside the Alfven surface is kept 
in solid rotation while flowing outwards along magnetic 
field lines. The plasma flow then becomes super-alfvenic 
and its angular momentum is conserved thereafter. The 
physical process responsible for maintaining the solid ro- 
tation while V < ca '^^ the magnetic tension, which 
works against the tendency for the field-lines to spiral 
backwards in the azimuthal direction (keeping them as 
straight as possible). Beyond the critical surface, the 
magnetic tension loses its efficiency. 

Of course, in the real sun the transition from rigid 
to non-rigid rotation is a continuous and smooth one. 
The azimuthal magnetic field does not suddenly change 
from straight to spiralled. But quantitatively, it all works 
out as if there was a sharp transition. Following this 
picture, one can define an average Alfven radius {ta) 
— an "effective lever arm length" . The specific angular 



momentum fiux rate is then 

l^n^{rl). (7) 

The resulting torq ue applied on the sun is 
(|Matt fc Pudritzl[2008h 

T = -Mflo(ri) . (8) 

The angular momentum per unit volume J of a parcel 
of solar wind plasma rotating with azimuthal velocity 
is 

Jw — prsinOv^ . (9) 

The angular momentum per unit volume crossing a 
surface element dA is then — J^VrdA. Integrating 
over a spherical surface of radius and assuming axi- 
symmetry translates into 

= 27rrQ / pVrV^ sin^ 9d9 . (10) 
Jo 

We then define the magnetic spin-down time-scale as 

Stsd = ^ , (11) 

where Jq is the sun's angular momentum. We estimated 

Stt f^'^^ 

Jq = — flo / p (r) r^dr « 1.84 x 10"*® g cm^ s"^ 

(12) 

(|Gilman et al.l 119891 : iStixl I2002D using a seismical ly cal- 
ibrated solar model fo r p (r) (a CESAM model, iMorell 
[19971 IBrun etall[200l) . 

The main difficulty now lies in the definition of (r^). 
Figure [9] shows the Alfven surface at two different in- 
stants of the cycle {t = 0,3 yr). This critical surface 
shows a regular shape for most of the activity cycle, be- 
ing close to spherical at most latitudes (especially, higher 
latitudes). At low latitudes, though, the Alfven surface 
approaches a more cylindrical shape. Some irregulari- 
ties appear as inward incursions as B vanishes in current 
sheets (there's a small but finite outflow, so v / ca ^ oo 
there). The most evident example of such an excursion 
corresponds to the equatorial streamer [e.g Figurel^l left 
panel). We defined here the average Alfven radius {ta) as 
the average cylindrical radii of the Alfven surface (that 
is, TA: the distance to the axis) weighted by the local 
mass flux sinSpv crossing the surface 

j r"^ sin 6pw ■rA(6)hde 
J sm 0\\pw\\dO 

The sections of the Alfven surface crossing heliospheric 
current sheets are discarded from the computation (that 
is, the inwards incursions described above). The average 
Alfven radius {ta) was found to be about 9 R© at the 
minimum and 2.2 Rq at the maximum (the time average 
being about 6.3 R©). Figure [10] (top panel) shows the 
average Alfven radius as a function of time during the 
whole solar cycle. The Alfven radius correlates well with 
the global magnetic field's amplitude for the most part 
of the cycle, but shows a negative correlation with the 
polar surface magnetic field (compare with Figure 2]) 

The values found for [ta) were then used to compute 
the solar spin-down time-scale (eq. [TTl with tq = {ta) in 
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eq. nop and the specific angular momentum flux (eq. [71 
using {r\) rather than (ta))- The results are displayed 
in Figure [10] (middle and bottom panels) as a function 
of time. The spin-down time-scale 5tsd due to the break- 
ing torque exerted by the solar wind was found to vary 
between 3.7 x 10^° yr and 2 x 10^^ yr (the time average 
being 4.7 x 10^-'^ yr). The quantities (r^i) and 6tsd show a 
strong negative correlation, meaning the "effective lever 
arm length" mostly determines the breaking efficiency. 
Note that the mass flux rate evolution throughout the 
cycle works against this trend ( cf. Figure [8]) : the highest 
mass outflow rate occurs at about the activity maximum, 
but its effect is overwhelmed by the large variations of 
(va) during the activity cycle. 

6. DISCUSSION 

We have shown that the coronal magnetic field's topo- 
logical variations due to the solar dynamo have a strong 
influence over the solar wind's velocity and associated 
mass and momentum fluxes. The ratio of open to closed 
magnetic flux at the sun's surface also varies. These ele- 
ments altogether define the variations of the Alfven radii 
(both in latitude and in time) and how effectively the so- 
lar wind carries angular momentum away from the sun. 

The average Alfven radius {va) was found to be about 
9 R© at the activity minima and 2.2 Kq at the maximum. 
The global poloidal magnetic field's amplitudes follow 
the same trend. Unexpectedly, though, the polar mag- 
netic field amplitudes are anti-correlated with the global 
field's amplitudes. The polar magnetic field (at the sur- 
face) happens to be stronger at the moments when the 
closed to open flux ratio is higher. As a consequence, the 
polar magnetic flux-tubes show higher expansion rates 
at these moments of the cycle (higher expansion factors 
imply stronger B-field decays). This behaviour could be 
related to the nature of the particular dynamo model 
used in this work, which is a simplified one (as exposed 
in i 32.1|) . We should nevertheless note that, in spite of its 
simplicity, the model reproduces well many solar coro- 
nal features. Namely, the topological characteristics of 
the poloidal magne tic field we obtain ed are consistent 
with those found bv lWang et al.l (|2003[ ) (using a different 
type of model for the surface field). Also, the latitudi- 
nal distribution of the wind velocities d uring the cycle we 
compu ted agree with those proposed bv lWang fc SheelevI 
(|2006D (well matche d by ULYSSES data; s ee their Fig- 
ure 3) and those by iTokumaru et al.l (|2010f) (built from 
IPS data collected during 2 solar cycles). 

The mass loss rate (as carried away by the wind) we 
computed averages out in time to ~ 5.3 x 10"^'* Mq/jt, 
a value comparable to in-situ measures near the Earth 
(giving a standard value of about lO"^** Mq/jt if the 
wind flow was spherically symmetric). Conversely, the 
solar spin-down time-scale due solely to the torque ex- 
erted by the magnetised wind outflow averages out to 
value of the order of 10^^ yr, above the standard 10^° yr. 

One reason for this disparity may lie in the fact that 
we used a one-fluid isothermal MHD wind model. Con- 
sidering a proper (i.e, non-isothermal) treatment of the 
energetics of the wind could lead to different quantita- 
tive estimat i ons o n the wind velocities. For example, 
iPinto et"al] ()2009[ ) found a broad range of wind veloci- 
ties in their configurations including a chromosphere - 
corona transition region and a coronal heating flux de- 



pendent on the local magnetic field's strength. The ac- 
tual asymptotic wind velocities depended strongly on the 
radial distribution of the mechanical heat fiux dissipa- 
tion. Under the right conditions, a non-isothermal two- 
fiuid app roach can also produce higher bulk velocities 
(see iEndeve fc Leeri[2001t IGrappin et al.| [2011'). An ad- 
ditional word of caution should be given in respect to 
the choice of a constant and uniform coronal tempera- 
ture. It is known that the mass flux in an isothermal 
wind flowing within a given magnetic flux-tube depends 
on the coronal temperature T as 

pv cx exp i-C/T) x / (r) , (14) 

Ccf. IParked 11961 ILeer fc HolzeJ [19791 IHansteen fc Leeil 
119951 aniong others) where po is the density at the base 
of the flux-tube, C a constant and f (r) oc A (r) /Ao a 
geometrical factor describing the radial expansion of the 
flux-tube. The absolute values found for the wind's mass 
flux are therefore very sensitive to the choice of the pa- 
rameter T, while the amplitude of its variations during 
the cycle (keeping T constant) are only due to the geo- 
metrical factor / (r) and the (small) variations in po. Ad- 
ditionally, it has been suggested that in non-isothermal 
scenarios the mass flux scales linear ly with the mechani- 
cal heat flux in j ected at the sur f ace (iLeer fc Holzeii[T979l : 
IHammei<[l982l: IWithbrod [I988I: IHansteen fc Leedll995D . 
We emphasize that, nevertheless, the results described 
in this manuscript depend essentially on other parame- 
ters of the problem, namely the cyclic variations in the 
global coronal topology and the consequent open flux- 
tube expansion factors (cf. the function / (r) above) . 
The isothermal approximation should in this case pro- 
vide a good assessment of the variability of the global 
wind properties during the activity cycle. Future work 
will address these issues by taking into account a more 
complete (i.e, non-isothermal) treatment of the energet- 
ics of the wind. 

Also, we neglected any effects due to the rotation on 
the wind's velocity (which should be small, anyway, as 
the sun is a slow rotator). 

The most remarkable feature found is the large tem- 
poral variability of these quantities (Alfven radius and 
spin-down time-scale) over an activity cycle. The sun 
loses mass and momentum non steadily, which translates 
into variable braking torques applied to the surface. We 
should add that the latitude and the width of the open- 
field and wind sources at the surface also vary in time, 
meaning the torque is applied at different places at differ- 
ent times of the cycle. During most of the activity cycle, 
the open-field crosses the solar surface at high latitudes 
only. The exception is the interval about the polarity in- 
version, when the polar coronal holes are greatly reduced 
and low latitude open-field fiux-tubes appear. We will 
consider whether this effect can be translated into mean- 
ingful upper boundary conditions (forcing, mass and mo- 
mentum sources/sinks) for the solar convection and dy- 
namo processes in future research. 

Now, as the solar cycle progresses the effects of an ever- 
changing magnetic field topology add up to the variations 
of the global amplitude of the field. This makes it hard to 
distinguish between their respective contributions. Let 
us first clarify the status of our coronal model in regard 
to this issue. We ran our model through 2 additional 
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Figure 11. Mean Alfven radius as a function of T. Top: variation 
of the amplitude of the external magnetic field for our ti = yr 
case. Middle: dipolar benchmark case (varations of the exter- 
nal field's amplitude, as above). Bottom: the same relation, but 
throughout the cycle. The labels ti to tg identify the instants 
shown in Figure \3\ The label ti appears both on the top and 
bottom panels, as it corresponds to the same state. 



test cases in which we varied only the external magnetic 
field's amplitude, and not its topology (i.e, the potential 
magnetic field only varied by a multiplicative factor be- 
tween consecutive runs). In the first case, we took the 
first element of our time-series (t = in the figures in 
the preceding sections) and varied the amplitude of the 
external field by a constant multiplicative factor. The 
second case corresponds to a purely dipolar external field 
B*^ whose amplitude we also varied, letting us perform 
a comparison with previous studies. We then computed 
the mean Alfven radius (r^) and the mass loss rate M 
for each new run and tested against the correlation 



M 



suggested bv iMatt fc Pudritd ()2008[ ). where ||B|| is the 
global field's amplitude. This scalar quantity is well de- 
fined for any given analytical external magnetic field but 
not for the ones obtained numerically from the dynamo 
model. We chose here to associate it to the unsigned 
magnetic flux integrated over the solar surface (after 
comparing both definitions for the dipolar test cases). 
Figure [TT] (top and middle panels) shows that the ex- 
pression above fits our 2 test cases for a power law in- 
dex m ~ 0.2, which is in good agreement with what 
IMatt fc Pudritz (2008) found for their dipolar case. This 
result validates the (rather large) variations in Alfven ra- 
dius we found. In addition, we plotted the same relation 
for the computations running through the whole cycle 
(Figure [11] bottom panel) . The same power law is still 
observed for most of the cycle, but it gets slightly steeper 
just before the polarity inversion and then deviates away 
from it (flipping sign) immediately after (between t = i 
and t — b yr), drawing a closed cycle in the diagram. 
This result suggests that the Alfven radius correlates 
positively with ||i?|p/M when the dipolar component of 
the open magnetic fleld is strong (which happens during 
most of the cycle), but may correlate negatively for some 
other higher order geometries. 

As expressed earlier in the text, we computed a suc- 
cession of steady-state wind solutions rather than the 
detailed dynamical evolution of a particular event. One 
justification for this choice is of numerical nature: the 
disparity of time-scales for the magneto-convective and 
coronal phenomena is very large. The second and per- 
haps more important justification is that varying the po- 
tential magnetic field in time introduces a non causal 
perturbation to the system. That would correspond to 
an instantaneous propagation of a perturbation to the 
field sources. This violation of causality remains even 
for a continuous and arbitrarily slow evolution of the 
background potential field. Note that localised pertur- 
bations to a MHD system are propagated away with 
phase velocities which correspond to the MHD wave 
modes. These can be arbitrarily small (even null in some 
places) and very anisotropic. To work around this is- 
sue one would need to include the field sources in the 
domain (and then self-consistently compute the system's 
response to their perturbation) . Alternatively, one could 
emulate the proper physical behaviour of such a system 
by injecting/propagating the hydro-magnetic perturba- 
tions through the numerical boundaries. This issue will 
be the subject of future work. 

7. SUMMARY AND CONCLUSIONS 

We have performed MHD numerical simulations cou- 
pling a solar dynamo model with a corona and solar wind 
model throughout a complete activity cycle. In short, we 
found that: 

• The latitudinal distribution of the asymptotic wind 
velocities is sensitive to the magnetic topology as 
it varies during the solar cycle. The fast wind - 
slow wind pattern s hows good qualit a tive a gree- 
ment with those in I Wang fc SheelevI (|2006f ) and 
iTokumaru eTall (|2010[ ). as shown in Figure 111 

• The polarity reversal happens rather abruptly in 
the corona, in contrast with the progressive evolu- 
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tion of the solar wind's velocities and of the surface 
magnetic field (Figure 14]). 

• Sun's global mass loss rate, Alfven radius and mo- 
mentum flux all vary considerably throughout the 
cycle (Figures m and [TU]) . The dominant causes are 
the position and latitudinal extent of the photo- 
spheric sources of solar wind and the geometry of 
the Alfven surface. 

• The zones of application of the braking torque due 
to the wind vary in time. Overall, the wind's break- 
ing torque should contribute to slow down the sur- 
face layers at high latitudes, but regions of appli- 
cation of torque appear occasionally at lower lati- 
tudes. 

Future work will focus on testing other types of solar 
dynamo models, using Babcock-Leighton flux transport 
based on large scale meridi onal circulation (jPikpati et al.l 
I2OOIJ J ouve &: Brun"2007') or turbulent magnetic pump- 
ing (jGuerrero fc dc Gouvcia Dal Pino 2008), as well as 
non-isothermal winds. We will furthermore consider the 
formation of a non-rigid coronal rotation and its effects 
on the solar wind properties. 
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